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Abstract. We are interested in the geometry of the group T> q (M) 
of diffcomorphisms preserving a contact form 9 on a manifold M. 
We define a Riemannian metric on T> q (M), compute the corre- 
sponding geodesic equation, and show that solutions exist for all 
time and depend smoothly on initial conditions. In certain special 
cases (such as on the 3-spherc), the geodesic equation is a simpli- 
fied version of the quasigeostrophic equation, so we obtain a new 
geodesic interpretation of this geophysical system. We also show 
that the genuine quasigeostrophic equation on S 2 can be obtained 
as an Euler- Arnold equation on a one-dimensional central exten- 
sion of Tid'Dq(M), and that our global existence result extends to 
this case. 

If E is the Reeb field of 9 and [i is the volume form, assumed 
compatible in the sense that div E = 0, we show that T> q (M) is 
a smooth submanifold of T>e,h{M), the space of diffeomorphisms 
preserving the vector field E and the volume form [i, in the sense of 
H s Sobolev completions. The latter manifold is related to symmet- 
ric motion of ideal fluids. We further prove that the corresponding 
geodesic equations and projections are C°° objects in the Sobolev 
topology. 



1. Introduction and Background 

Contact geometry is the odd-dimensional analogue of symplectic ge- 
ometry jEKMl [G], and the purpose of this article is to extend some 
results on the Riemannian geometry of symplectic forms |E2] to the 
contact case. On a manifold M of dimension 2n + 1, a contact struc- 
ture is a distribution of 2n-planes which maximally does not satisfy the 
Frobenius integrability condition at any point. If the contact structure 
is orientable (as in typical examples), it is determined as the null space 
of a contact form 8, a 1-form on M for which maximal nonintegrability 
translates into 
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everywhere. On M 2n+1 with coordinates (x 1 , . . . , x n , y 1 , . . . , y n , z), the 
basic example is given by 



By the Darboux theorem [G], we can always find local coordinates on 
any manifold M such that 9 is given by (j2J). 

It is important to note that if 9 is a contact form, then for any 
nowhere zero function F on M, the 1-form F9 also satisfies ([1]) and 
its null space is the same contact structure. Hence if we care primar- 
ily about the contact structure, the symmetries are the group V g (M) 
of diffeomorphisms t] such that t]*9 = F9 for some positive F; if we 
care primarily about the contact form, the symmetries are the group 
T> q (M) of diffeomorphisms r] with t]*9 = 9. In this article we will 
work with the latter group, which following Ratiu-Schmid |RSj . we call 
the quantomorphism group. In a separate article [EbPrj we will study 
the geometry of the contactomorphism group T>g(M), which also has 
a number of interesting features and has not been studied from the 
Riemannian point of view at all. 

The quantomorphism group was studied by Omori [O] (as an ILH 
group), Ratiu-Schmid [RS] (as a Hilbert manifold), and Smolentsev [S2] 
(as a formal Riemannian manifold). More recently several authors have 
observed that the quantomorphism group is the proper configuration 
space for the geodesic Vlasov equations (our geodesic equation is es- 
sentially a special case of this). Holm-Tronci |HT] proposed a general 
quadratic Hamiltonian on the quantomorphism group and studied its 
moments. Gay-Balmaz and Vizman |GVj computed the momentum 
maps for this equation, and Gay-Balmaz and Tronci [GTj found an 
interesting totally geodesic subgroup. We would like to thank Darryl 
Holm, Tudor Ratiu, and Frangois Gay-Balmaz for very helpful discus- 
sions. Our main contribution here is local and global well-posedness of 
the geodesic equation in Sobolev spaces. 

We now review some well-known properties of contact structures. 

Proposition 1.1. There is a Reeb vector field E which is uniquely 
defined in each tangent space by the conditions 9(E) = 1 and tE,d9 = 
0. The linear map u y i u d9 from T X M to T*M yields a pointwise 
isomorphism 7 from the null space of 9 to the annihilator of E. 

Proof. In Darboux coordinates we have d9 = J2k=i A dy k , and it 
is easy to check that the Reeb field must be E = d z . The null space 
of 9 is spanned by the basis {d x k,d y k — x k d z : 1 < k < n}, and 
we have 7(<9 x fc) = dy k and ^(d y k — x k d z ) = —dx k for all k. Since 
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dy k (E) = dx k (E) = 0, the map 7 takes a basis of vectors annihilated 
by 9 to a basis of covectors that annihilate E. Since these definitions 
are coordinate- invariant, they define global objects E and 7. □ 

Properties of the Reeb field determine the structure of the quanto- 
morphism group, in the following way. 

Proposition 1.2. Any diffeomorphism rj for which rj*9 = 9 must also 
satisfy rj*E = E. As a result the quantomorphism group T> q (M) is a 
subgroup ofT>E{M), the group of diffeomorphisms preserving E (alter- 
natively, the diffeomorphisms commuting with the flow of E). 

Proof. We want to show that rj*9 = 9 implies that Dt] o E = E o rj. We 
write r]*E = Dt] o E o rj~ l so that rj*E is another vector field on M. 
Then we have 

L v ,Ed9 = L E {rfd9) o r/ _1 = i E (d9) o rj' 1 = 0, 

along with 

9{r)*E) = rf9{E) of^lo rj' 1 = 1. 
Since 77* E and E satisfy the same conditions, they must be equal. □ 

If \i is the Riemannian volume form, and we assume (as we shall 
always do) that fi is a constant multiple of 9 A (d9) n , then clearly 
ri*9 = 9 implies that rfii = \i. Hence Proposition 11.21 shows that 
Vg(M) is a subgroup of V E41 (M), which we refer to as the group of 
volumorphisms with symmetry. 

Example 1.3. Our most important example will be S 3 , considered as 
a group of unit quaternions in IR 4 . It has a global basis of left-invariant 
vector fields {E\, E2, E%} satisfying the relations 

[Ex, E 2 ] = -2E 3 , [E 2 , E 3 ] = -2Ex, and [E 3 , E x ] = -2E 2 . 

Letting {a; 1 ,^; 2 ,^; 3 } denote the dual basis, a contact form is given by 
9 = to 1 , with Reeb field E = E\. We have d9 = 2cu 2 A to 3 , so that 
9 A dO is twice the Riemannian volume form. If we declare {E{\ to be 
orthongal with ||-Ea|| = a and H-E2II = H-E3II = 1 for some parameter 
a > 0, then we have the Berger metric on S* 3 . (See for example [F].) 

The Reeb field E — E\ is Killing, and all orbits are closed and have 
the same length. The quotient is S* 2 , and the projection n: S 3 — > S 2 
is the well-known Hopf fibration. The pullback of the standard area 
form on S* 2 is a constant multiple of d9. Any quantomorphism 77 on 
1S 3 generates a symplectomorphism £ on S 2 by the formula ((p) = 
7r(?7(7r~ 1 (p))), which is well-defined since 77 commutes with the flow of 
Ei and hence is constant on the fibers of n. Notice that the flow of E\ 
gives a family of quantomorphisms which preserve the fibers of ir, and 
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hence all map to the identity symplectomorphism of S 2 . This is the 
typical behavior we expect quantomorphisms to have; see below. 

Not every contact manifold has an interesting quantomorphism group. 

Example 1.4. On M = T 3 , we can check that 9 = sin zdx + cos zdy 
satisfies 9 A d9 = dx A dy A dz, so that 9 is a contact form. The Reeb 
field is E = sin z d x + coszd y . Every quantomorphism must preserve 
the Reeb field, but the Reeb field has nonclosed orbits whenever tan z 
is irrational, and hence any function which is constant on the orbits 
must actually be a function only of z. It is then easy to see that the 
identity component of V q (T 3 ) consists of diffeomorphisms of the form 

rj(x, y, z) — (x + p(z) sin z + p'{z) cos z, y + p(z) cos z — p'{z) sin z, z) 

for some function p: S 1 — > R. This group is abelian, so any right- 
invariant metric will actually be bi-invariant, and all geodesies will be 
one-parameter subgroups. 

Instead what we want is for all the orbits of the Reeb field E to 
be closed and of the same period, so that the flow t i— > <p(t, x) of E is 
periodic. In this case we say that the vector field and the corresponding 
contact form are regular [RS] . When this happens, there is a Boothby- 
Wang fibration |BWj n: M — > N, where N is a symplectic manifold 
of dimension 2n with symplectic form u satisfying tt*u = d9. Using 
Darboux coordinates (x k , y h , z) on M, we get coordinates (x k , y k ) on 
iV such that ui = Y^k dx k A dy k . The Hopf fibration of Example 11.31 is 
the basic example. 

Example 1.5. To get a Boothby-Wang fibration, it is important that 
the orbits are not only closed but of the same length. Here we give 
an example of what happens when this fails. The Hopf fibration is 
commonly constructed by considering S 3 as the unit sphere in C 2 . Using 
complex coordinates (w, z) we let the circle action be r i— > (e lT w, e tT z) 
so that all orbits are closed and of period 2tt. A variation of this, which 
also comes from a contact form, is the action r i— > (e tkT w, e llT z) where 
k and I are relatively prime. In this case the orbits have period 2n 
except when w or z is zero in which case they have period of 2n/k or 
2ix/l. Thus the quotient 5 3 /5' 1 has singularities and is not a manifold 
at the images of these points. This construction actually comes from 
a simple problem in mechanics. Take two harmonic oscillators with 
periods 2ir/k and 2n/l respectively. Each of their motions is a curve in 
C so together they are in C 2 . S* 3 C C 2 will then be a constant energy 
surface and the motion on it will be the S 1 action described above. We 
thank H. Hofer (personal communication) for this example. 
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Let us now review some facts about the group of quantomorphisms. 
At first we will assume that all objects are C°°; later we will work 
with spaces of Sobolev diffeomorphisms, so that we can do analysis 
on Hilbert manifolds and use tools of differential analysis such as the 
inplicit function theorem. 



Proposition 1.6. The tangent space to the quantomorphism group 
V q (M) = {rj E V(M) | ri*9 = 9} at the identity is 



For any f with E(f) = 0, the field u = Sgf satisfies C u 9 = 0; con- 
versely if C u 9 = 0, then u = Sgf for some f with E(f) = 0. 

Proof. Since rj G T> q (M) if and only if 7]* 9 = 9, a tangent vector u 
at the identity must satisfy C u 9 = 0, or d{9{u)) + i u d9 = 0. Letting 
/ = 9(u), we have i u d9 = —df. Applying both sides to the Reeb field 
E, we have E(f) = d9(E, u) = 0. Since E(f) = 0, we know df is in the 
annihilator of E, and since 7 is an isomorphism from the null space of 9 
to the annihilator of E, we get a unique field v for which 9(v) =0 and 
7(1?) = —df. Then we must have u — v + fE in order to satisfy both 
9(u) = f and j(u) = —df. The rest of the proof is straightforward. □ 

Note that any C°° function /: M — ¥ M satisfying E(f) = is con- 
stant on orbits of E and thus may be viewed as a function / : iV — > K 
where N is the symplectic quotient. We can thus identify the space 
jFe{M, R) from (j3J) with J r (A r , R) whenever convenient. Also since u 
is a symplectic form on N we get an isomorphism to' 9 : TN — > T*N 
defined by cu^(v) = l v uj. We denote its inverse by uA 

Using Sg we can define a type of Poisson bracket on Te{M, R). We 
define 



(3) 



TuD q = {S e f\fe T E (M,R)}, 



where 



(4) T E {M, R) = {f:M^R\ E(f) = 0} 




(6) 



{f,g} :=Sgf\g). 



Thus in Darboux coordinates we have 
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and 



II 



df dg df dg 



(8) 



dx k dy k dy k dx k 



k=l 



(since of course |f = for g e T E {M, M)). 

Proposition 1.7. Sg : J- E (M, M) — > T i( fP q is an algebra isomorphism; 
that is, it takes Poisson brackets into Lie brackets. 

Proof. Let u = Sgf and v = Sog. Then since Sg is surjective, there 
exists h e J-£;(M, M) such that Soh = [u,v]. Then 9(Sgh) = h and 



9([Sgf,S e g]) = -d9(Sgf,Sgg) + Sgf(9{S e g))-Sgg(9(Sgf)) 

= -L Sef d9(S e g) + S g f(g) - S e g(f) 

= df(S e g) + {f,g}-df(S e g) 

= U,g} 



Note that ir*(Sgf) = —oj$(df), and the identification of ^(M, R) 
with F^NjM.) takes the Poisson bracket of J~e(M, M) into the Poisson 
bracket of F(N, R). Also since the standard bracket has the property 



as well, where [i is the contact volume form. This fact will be needed 
for the proof of Theorem 14.11 

If we denote the operator u i— > C u 9 by Lg, then by construction we 
have LgoSg = 0. We can then view the operators So and Lg as forming 
a short exact sequence if we choose their domains and ranges correctly. 
This will be important when we discuss smoothness of the geodesic 
equation, so we record the formal result. 

Proposition 1.8. Let T^iMi ^) denote the space of smooth functions 
f-.M^R such that E(f) = as in (0}. Let TgA^^M) denote 
the space of one- forms (3 on M such that le(3 = 0, ig<i/3 = and 
9 A d(3 A (d9) n ~ 1 = 0. Finally, let T id V E ^(M) denote the tangent space 
of V 'e ,h{M) , which consists of vector fields u such that divu = and 
[E,u] = 0. If Lgu = C u 9 and Sg is defined as in then 



Thus h = {f,g} so S g ({f,g}) = [S e f,S e g]. 



□ 




where v is the symplectic volume form on M, we get 




(9) 



-)• F E (M,R) ^ T ld V E ,^M) H T 9 A E JM) -»• 
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is a short exact sequence, i.e., the image of every map is the null space 
of the next map. 

Proof. That Sg is one-to-one from F E (M, R) is easy to check. The 
fact that Sg maps into T id V EfJi follows from Proposition II. 2[ as well as 
the fact that the volume form /i is a multiple of 9 A (d6) n . To show 
exactness at T id V Etfl (M) we pick u G T id D E ^(M). L e u = if and only 
if i u d6 = —d6{u). Letting / = 9(u), we find that u = Sgf so exactness 
here follows. 

Finally the fact that L e maps T id D E ^(M) onto T e A E41 (M) follows 
from the following computation. Let /3 G TgA E ^(M), and define u to 
be the unique vector field satisfying 6{u) = and j(u) = (3. Then 
Lg{u) = (3, and we just need to check that u G T- ld D E ^(M). We have 

£[e,u\9 — C E C U 9 — £>u£>e® = £eP = 0, 

and since 9([E,u]) = E{6{u)) - u(6(E)) - d6(E,u) = 0, we conclude 
that j([E, u]) = as well. Since the only vector field v satisfying 
6{v) = and j(v) = is zero, we conclude [E, u] = 0. The fact that u 
is divergence-free follows from the computation 

C u (6 A (d6) n ) = CJ A {d9) n + n6A dCJ A (d0) n ~ l = 

since C u 6 = (3. □ 

2. SOBOLEV MANIFOLD STRUCTURES 

In this section we discuss the Riemannian metric on T> q and its ge- 
odesic equation. We will show that T> q is a smooth submanifold of 
T> E pjl and that the geodesic equation is a smooth ordinary differential 
equation on T> q . The results of this section generally do not depend on 
the choice of Riemannian metric on M; we can always arrange things 
so that there is some metric on M which makes E a Killing field, by 
the regularity assumption on E. 

To discuss smoothness, we extend our spaces of C°° maps to Sobolev 
H s maps, where s is an integer strictly larger than dimM/2 + 1 = n + | 
(in order to ensure that any such map is C 1 ). These structures have pre- 
viously been studied in this context by Omori [D], Ratiu-Schmid |RSj . 
and Smolentsev [S2j . Our interest is in the quantomorphism group 
T> s q (M) as a submanifold of the volumorphism group with symmetry, 
T> S E AM), which in turn is a submanifold of T> S ^(M) (which is known to 
be a smooth submanifold of V s (M) by |EM| ). 

We will be dealing with many Hilbert manifolds, all of which are 
topological subgroups of the full diffeomorphism group V s (M), and we 
would like to show they are all submanifolds of this group. This is 
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not automatic even for nice diffeomorphism groups; for example the 
contactomorphism group 2?|(M) is a C°° Hilbert manifold but not a 
C°° Hilbert sub manifold of V s (M) (see Omori [O]). We will generally 
prove our subgroups are submanifolds either by constructing an explicit 
coordinate chart or by using the implicit function theorem. In order to 
avoid duplication, we will do this only in the simplest cases, and derive 
the other cases from it. 

Lemma 2.1. Suppose L and M are both C°° Hilbert submanifolds of 
a C°° Hilbert manifold N, and that L is a subset of M. Then L is a 
C°° Hilbert submanifold of M. 

Proof. For general properties of Hilbert manifolds, see [AMR] or [Lj. 
We will use the fact that a subset is locally a smooth submanifold if 
and only if the inclusion map is a smooth immersion. Let i : L — > M, 
j : M — > N, and k : L — > N be the inclusions. We know j and k are 
smooth immersions and want to prove that i is as well. Clearly joi = k, 
and if we knew that i were smooth, the Chain Rule would imply that 
it would have to be an immersion since Tj o Ti = Tk, with Tj and Tk 
both injective. 

Smoothness of i follows from the "universal mapping" property of 
submanifolds; see Lang jE]. More explicitly, we get smoothness using 
coordinate charts. Since M is a smooth submanifold of N, for every 
p G L C M there is a chart tp: W — >■ G from a neighborhood W of p 
in N to a Hilbert space G, such that for some closed subspace F C G 
we have 

i[)[WnM] = ip [W] fl F x {0}, 
where {0} denotes the zero element of the orthogonal complement F- 1 . 
The map -0 = irF°il>\wnM '■ WC\M — > F given by restricting the domain 
and projecting the range onto F is a coordinate chart on M. 

Now since ip is a coordinate chart on N, we know ip o k : L — > G is 
C°°. Thus Tip o if; o k: L — > F is also C°°. But 7Tpo?/;o&; = ?/>oz, and 
since ipoi is smooth for a coordinate chart i/), the inclusion map i must 
be smooth by definition. □ 

To begin we consider the set 

(10) V E {M) = {rjE V S {M) | 7)*E = E} 

of H s diffeomorphisms of M which preserve the Reeb field (or equiva- 
lently, the set of H s diffeomorphisms which commute with the flow of 
the Reeb field). We will prove this is a submanifold of T> S (M), a special 
case of a result due originally to Omori [O]; we provide the proof to 
make the paper somewhat more self-contained. 
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Theorem 2.2. Let E be a vector field on M with closed orbits all of 
the same period. The set T> S E (M) of H s diffeomorphisms preserving E 
(or commuting with its flow) is a C°° submanifold ofT> s (M). 

Proof. We may assume without loss of generality that E is a Killing 
field of some Riemannian metric on M; we can simply take an arbitrary 
metric and average it over a period of the flow of E. 

The coordinate charts on V s (M) are defined in terms of the Rie- 
mannian metric on M as follows: for any 77 G T> s we consider the linear 
space H*(TM) of H s vector fields V over 77, i.e., such that V(x) G 
T v i x \M for every x G M. Then using the Riemannian exponential map 
exp on M, we define an exponential map Exp^: H*(TM) — > V s by 
Exp^fV) = x H- exp v ( x j(V(x)) (i.e., for each x we follow the geodesic 
starting at position rj(x) with velocity V(x) for time one). If the H s 
norm of V is sufficiently small, then by the Sobolev embedding theorem 
the C° norm of V will be smaller than the injectivity radius of M, and 
thus Exp, is invertible on a small neighborhood Q of zero in H^(TM); 
its inverse is the desired chart. 

To prove D S E is a submanifold, it is sufficient to show that in these 
charts, for any 77 G V S E (M), we have 

(11) Exp, [nn if'] =Ex Vri [n}nv s E 

for some closed subspace of H*(TM). Set 

K s v = {vo V \[E,v] = 0} = R v [K^ 

where R v is the right-translation. Since K? d is the null-space of the 
continuous map v 1— y [E, v] from H s vector fields to if s_1 vector fields, 
it is a closed subspace of H? A (TM); since the right translation R v is an 
isomorphism, each space K!L is closed in H^(TM). 

Let ( T denote the flow of E; since E is a Killing field, every ( T is an 
isometry. Since isometries preserve geodesies we have, for any 77 G T> s 
and V G H*(TM), that 

(12) exp Cr(j?(:c)) ((( T )*V(x)) = Cr(exp„ (x) (V(a;))) 

for all x G M and r G M. Now if G then V = v o 77 for some 
t> G J^f d , and since [E, v] = 0, the flow £ r of E preserves v, i.e., we have 
{Ct)*v(x) = v(( T (x)) for all x G M and all r el, We conclude that 
(Ct)*V(x) = v (Cr(f]( x ))) ■ Combining with (fP2l) we obtain 

(13) eX PCr(r;0r)) { V {CAV(X)))) = Cr ( eX Pr?(;c) (^(t/(x)))) 

for any 77 G and f G i^f d . 
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Since 77 G V S E , we know that 77 o ( T = ( T o 77, so that ( TT3l) becomes 

Exp^V) o C T = Exp^V o C r ) = Cr o Exp^V) 

for any r G M. Since o := Exp v (V) commutes with the flow ( T of 
E, differentiating with respect to r at r = gives a*E = E, and we 
conclude that Exp^ actually maps fl fl into V S E (which proves that 
the left side of (fTTT) is a subset of the right side) . 

To prove the opposite inclusion, let a be any diffeomorphism in 
Exp r) [fi]nX>| ; ; we want to show that a = Exp^fV) for some V G Oflif*. 
We already know that Exp^ is a diffeomorphism on Q, so that there is 
a unique V G fl with Exp v (V) = a; thus we just need to show that 
V G K^, i.e., that V = v o 77 for some t> with [E, v] = 0. 

Since a G we know that cr o ( T = ( T o <x; hence for each x G M 
we have 

( T (exp v{x) (V(x))) = exp v{CAx)) (V(Cr{x))) = eW( T ( v (x)) {V{Ct{x))), 

using the fact that 77 G V S E as well. Since isometries preserve geodesies, 
we also know that 

Cr(ex Pr?(x) (y(x))) = exp Cr(r?(:r)) ((C r )*F(x)). 

Now by our construction of f2, we know the supremum norm of V is 
smaller than the injectivity radius of M, so that the minimizing geo- 
desic between r](x) and o~(x) is unique for each x. Since both V(( T (x)) 
and (Ct-^V^x) are initial velocity vectors for the geodesic at ( T (v( x )) 
which at time one reaches Cr(7(^)), they must be equal for all suf- 
ficiently small t. Since V — v o rj, we conclude that v ^ T (rj(x))) = 

(Ct)* v (v( x )) f° r & H x i i- e -; that f (Ct(?/)) = (Ct)*^(?/) f° r ah 1/ £ M. Now 
differentiating with respect to r at r = 0, we conclude that [E, v] = 
as desired. □ 

We now consider the set of H s diffeomorphisms T> S E M (M) preserving 
both the Reeb field E and some volume element \x. For this to be in- 
teresting, we want the Reeb field itself to preserve the volume element, 
so we assume that E is divergence- free. Of course this will be the case 
if E is Killing. We will prove that T> E JM) is a C°° submanifold of 
T> E (M); this is proved in Smolentsev |Slj using a different technique. 

Theorem 2.3. Let M be a compact manifold. Suppose E is a smooth 
vector field with closed orbits all of the same period, and let \i be a 
volume form invariant under the flow of E (so that E is divergence- 
free). Let T> S E n(M) denote the set of H s diffeomorphisms rj such that 
rj*E = E and rf [i = \i. Then T> E is a smooth submanifold ofD s E . 
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Proof. The easiest way to proceed is to use the implicit function theo- 
rem. Consider the function V given by V(rj) = r\*\i- If 77 ET> S then 77*// 
is an i? s_1 n-form, given explicitly in coordinates by the Jacobian de- 
terminant of 77; since this just involves multiplying the first derivatives 
of 77 together, and since multiplication of H s ~ x functions is smooth jE], 
we know that V is a C°° function from V s (M) to H s ~\A n ). 

However we want to consider V on the P-invariant diffeomorphisms, 
which will necessarily map into the P-invariant n-forms. Since T> S E is a 
smooth submanifold of V s , the restriction of V is smooth on D S E . If ( T 
again denotes the flow of E, our assumptions imply that = /1, so 
that if 77 G T> S E and fi = r]*fi, we must have 

= (Cr)*(r/» = V*(Cr)*fi = r)*n = L~L. 

We can write any such H s ~ x volume form fx as fx — Ffi for some iP _1 
function F such that J M F/x = f M fx- Also since fx is ^-invariant, the 
form fi will be £"-invariant if and only if the function F is. Let 

H E - 1 \A n ) = \fteH E - 1 [ fi = [ ,A, 
I Jm Jm ) 

a space of E- invariant n-forms. We can identify this with H E ~^(M, M), 
the space of i?-invariant functions with mean one. (The fact that 
any such function is the image of some diffeomorphism rj follows from 
Moser's result |Moj . as generalized to the Sobolev case in |EMj .) One 
can prove exactly as in Theorem O that H S E ~\M, R) is a C°° sub- 
manifold of H S ~ 1 (M, R), the space of all .£P -1 functions on M; thus in 
particular it is a C°° Hilbert manifold. 

We thus have a C°° function V from the Hilbert manifold T> S E to 
the Hilbert manifold H E ~^(A n ), and the inverse image V~ 1 {li} is the 
set T> E[i \ thus to prove T> E[i is a submanifold it is sufficient to show 
that the derivative of V is surjective in each tangent space. We have 
{TV\{v o 77) = r]*(C v Li) for any 77 G V S E and v G Ti,fD E , and since 77* 
is an isomorphism of H s ~ l spaces of n-forms, it is sufficient to show 
that (TV)id is surjective as a map from T id V s E to T IM H E ~ 1 (A n ). Since 
every fi G H E ~^(A n ) has the same integral as /x, this tangent space 
consists of the P-invariant _£P -1 n-forms which integrate to zero over 
M. This space can be identified with H e ~q(M, R), the space of E- 
invariant functions whose integral is zero. 

So let A be an H s ~ 1 P-invariant n-form with f M A = 0. Then A = ffi 
for some / G H e ~q((M, R). Choose a Riemannian metric on M such 
that /i is its Riemannian volume form, and integrate over the orbits of 
E to obtain a metric such that E is a Killing field; since E preserves 
fi, the new metric will still have fi as its volume form. Using the 
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Hodge decomposition, we can write / = divt>, where v = VA -1 / is 
an H s vector field. A direct calculation shows that C v fi = A. Since 
E is a Killing field, it commutes with the Laplacian, and we have 
[E,v] = V£(A-7) = VA^E(f) = 0, so that v G T id V s E as desired. 
Thus (TP) id is surjective, and the implicit function theorem implies 
that T> E is a smooth submanifold of T> S E . □ 

Corollary 2.4. Under the assumptions of Theorem \2.3[ the space 
V E (M) is a C°° submanifold ofV^iM), the space of H s diffeomor- 
phisms preserving the volume form fi. 

Proof. Combining Theorem 12.31 and Theorem 12.21 we see that D S E is 
a smooth submanifold of V s . Since T>% ,, is a subset of T>f,, and T> s „ is 
a smooth submanifold of V s by [EM] . Lemma [2.11 gives the result. □ 

This Corollary will be more useful to us than Theorem I2.3[ since it 
allows us to obtain a connection on T> S E directly from the one on T> s 
constructed in |EM] rather than having to construct it by imitating the 
proof there. 

Now we discuss the quantomorphism subgroup. Ratiu and Schmid 
[RS] proved that this is a C°° submanifold of T> E (M) (the space of H s 
diffeomorphisms preserving the Reeb field E); we will use a modified 
technique to prove it is also a C°° submanifold of T> E AM). The latter 
space is more useful geometrically since it has a natural right-invariant 
L 2 Riemannian metric, while the natural L 2 metric on T> S E (M) is not 
right-invariant. 

The following Lemma is useful in the proof. 

Lemma 2.5. Suppose E is a Killing field on a compact Riemannian 
manifold M (without boundary) with closed orbits. Let (3 be a 1-form 
with Hodge decomposition 



where <fi is a function and u is a 1-form with 5u = 0. If (3(E) = and 
L E df3 = 0, then the same is true for each term: d<p(E) = 0, u(E) = 0, 
and L E du = 0. 

Proof. By Cartan's formula, C E a = L E da+d(i E a) for any 1-form a, and 
hence Ce commutes with d for any vector field E. Typically Ce does 
not commute with 5; however if E is Killing it does, because 6 is the 
adjoint of d. Hence for any k- forms (3 and a we have J M g(CE(3, o~) d\x = 
— J M g((3, Ce<j) dfi since Ce9 = for the induced metric on Axforms 
and CeH = for the Riemannian volume form. Hence Ce is skew- 
self adjoint, so that 



(14) 



(3 = d(j) + uj 



6c E = -d\c E y 



(dc E y 



C* E d* = C E 5. 
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By assumption we have CeP = 0, and so we have = d(CE(p) + ^E<^', 
since SCeoj = Ce{5oj) = 0, this must be the Hodge decomposition of 
the zero 1-form. But the Hodge decomposition of a 1-form is unique 
so dC E (p = and C E <^ = 0. Thus £e</> must be a constant. But 
Ce4> = E(<j>), and since E has closed orbits, the only way E(<f>) can 
be constant is if E(<p) = 0. Now applying both sides of (Tl4|) to E we 
obtain uj(E) = 0. Finally since ^gu = and u(E) = 0, we must have 
L E du = 0. □ 

Now we prove T> s q (M) is a smooth submanifold of V S E (M). Ratiu 
and Schmid [RS] gave a proof of this with a minor error, which we fix 
using Lemma [2.51 

Theorem 2.6. Let M be a contact manifold of dimension 2n+ 1 with 
contact form 9 and regular Reeb field E, and assume the volume form 
fi on M is a constant multiple of 8 A (d6) n . Then for s > n + 3/2 7 the 
space T> q (M) of H s diffeomorphisms preserving 9 is a C°° submanifold 
of the space T> S E (M) of H s diffeomorphisms with symmetry. 

Proof. We use the implicit function theorem for Hilbert manifolds. By 
definition the quantomorphism group V q (M) is the set of r] such that 
71*6 = 6, so it is tempting to simply use the function rj i— > rj*6 in the 
implicit function theorem; however there are a couple of subtleties that 
arise, primarily due to the fact that although rfQ is only in H s ~ x for 
an H s diffeomorphism t], we also have drj*8 = rj*d9 which is in H s ~ l . 
Thus 7]*9 is smoother in some directions than in others, and we need 
to keep track of this in order to show that 7] h- 7]*9 is a submersion. 

Let A S E ~ 1 (M) denote the set of all ordered pairs (a, da) satisfying the 
following four conditions: 

(15) aeiF-^A 1 ) and da G H S ^ 1 (A 2 ); 

(16) a(E) = 1 and lEda = 0. 

These conditions are motivated by the fact that they are satisfied for 
any 1-form rfQ for rj £ T> S E „(M), as we will show in a moment. First 
we will show that A^T 1 is a smooth Hilbert manifold. Consider the 
space if s ~ 1 (A 1 ) x H S ~ 1 (A 2 ) (the space of ordered pairs of H^ 1 1- 
forms and H s ~ l 2-forms on M, not necessarily related to each other). 
The subset {(a, da) \a £ H s ~ 1 ,da £ H 13 " 1 } given by condition (1T51) 
defines a closed linear subspace, and since the maps a —> a(E) and 
da — > LEda are linear and continuous, condition ffl6|) defines a closed 
affine subspace of that subspace. Thus A^T^M) is a C°° affine manifold 
with tangent space 

T a A s -\M) = {(p, dp) £ H'-^A 1 ) x H S ~\A 2 ) | P{E) = 0, i E dp = 0}. 
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Now let F : V S E (M) -> if s_1 (A 1 )xif s_1 (A 2 ) denote the map F e {rj) = 
(rj*9,rj*d9). In coordinates such that 9 = J2k ^k( x ) dx k , the 1-form of 
Fgijf) is given by 

(r ] *9)(x) = J2^(v(*))^dx>- 

thus if rj is H s , this is an H 8-1 1-form; similarly d(rj*9) = rfd9 is an 
jjs-i 2-form. Since multiplication and composition with C°° functions 
are smooth for H s ~ x functions with s — 1 > dim M/2 [E], Fa is a C 00 
map. Since any r\ G V S E (M) preserves the Reeb field, it is easy to see 
that (p = r/*9 satisfies <f>(E) = 9(E) = 1 and t^o^ = iEd9 = 0. Thus 
the image of Fg is contained in the space A^ . We will show that Fg 
is locally surjective onto it. 

By definition of the Lie derivative, the derivative of Fg is the map 
(TFg) v (u o rj) = r]*C u 9 for any V G T> E JM) and any u G T id V s E ^(M). 
Thus to show that TFg is surjective everywhere it is sufficient to show 
it is surjective at the identity, since 77* and composition with r\ are both 
isomorphisms. 

Let (J3,dp) G TgA s E ~ 1 (M); we want to find u G T id V s E (M) such that 
C u 9 = (3 and therefore C u d9 = d/3. In other words we want u = fE + v 
where / G H S (M, R),ve H S (TM), E(f) = 0, 9(v) = and [E, v] = 0, 
such that (3 = C u 9 = df + i v d9. The fact that we have to differentiate 
/ and not v in this formula is what makes this a bit tricky. Since /3 is 
in if s_1 , we can write /3 = d<j) + u where is a function in H s and uj 
is an i/ s_1 1-form with 5u = 0. However since we know d/3 is also in 
H s ~ l , and since we have d/3 = du = (d + 5)u, we get <5c?/5 = Au. We 
conclude that is actually in H s . Also by Lemma 12.51 we conclude 
that 

E{4>) = 0, u(E) = 0, and L E du = 0. 

The first formula implies that we can choose / = </>. The second formula 
implies that there is a unique vector field v such that 9{v) = and 
L v d9 = u, because the map 7 defined in Proposition II .11 is invertible as 
a map from the null space of 9 to the annihilator of E at every point. 
Since d9 is smooth, we see that v is also in H s since u is. 

Finally, the formula Lj^du = implies that [E,v] = 0; this comes 
from combining the fact that 9([E, v}) = E(9(v))-v(9(E))-d9(E,v) = 
with the fact that C[e, v ]9 = CeC v 9 — C v Ce9 = Ceoj = to obtain 
L [E,v]d9 = 0. Since [E,v] is in the null space of both 9 and d9 at every 
point, it must be zero everywhere. 
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We have thus shown that Fq is a submersion from the smooth mani- 
fold V S E {M) to the affine space A^ _1 (M), and thus V s q (M) = F" 1 ^, dO) 
is a smooth submanifold of T> S E (M). □ 

Corollary 2.7. Under the same assumptions as Theorem \2.6l the 
quantomorphism manifold T> s q (M) is a C°° submanifold ofT> s Etl (M), 
the space of volume-preserving maps preserving E. 

Proof. This follows from Lemma [2. II □ 

3. Connections and geodesic equations 

In the last section we were concerned only with the different iable 
structure of the diffeomorphism groups. Having established smoothness 
of all the basic objects, we now investigate the (weak) Riemannian 
geometry. First we define weak metrics, then we establish the existence 
of a smooth Levi-Civita connection, and finally we derive the geodesic 
equation. Smoothness of the connection guarantees that the geodesic 
equation is a smooth ODE on a Hilbert manifold, and hence we have 
a smooth Riemannian exponential map (and in particular, local well- 
posedness of the geodesic equation). 

We start with the standard L 2 metric on T> s (M), defined for vectors 
U,V e T V V S by 

(17) ((U,V)) V = [ (U(x),V(x)) v(x) dfj,(x) 

in terms of a Riemannian metric on M and its Riemannian volume 
form /i. Although this metric is neither right-invariant under the group 
action nor strong enough to generate the Sobolev topology, it is conve- 
nient and natural since its geodesies are given in terms of the pointwise 
geodesies on M. That is, for any rj ET> S and U G T V T> S , the geodesic 
on V s through rj with initial velocity U is 

(18) Exp v {tU) =x^ exp v(x) (tU(x)). 

The corresponding Levi-Civita connection on TT> S is C°°, as shown in 
[EM] . For right-invariant vector fields U and V on V s (M) such that 
U v = uorj and V v = v o rj for vector fields u and v on M, the covariant 
derivative is given by the expected formula Vy V = (V u v) o r\ where V 
is the covariant derivative for V s and V is the covariant derivative of 
M. Given the connection on V s , we will construct connections on sub- 
manifolds using orthogonal projections, and if the projection is smooth 
then the resulting connection will be smooth on the submanifold. 

On the smooth submanifold T> s ^ of volumorphisms, the restriction of 
the Riemannian metric (fTTI) is right-invariant under the group action. 
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The tangent spaces are given by 

TrPl = { U or]\dwu = 0}, 

and right-invariance implies that the orthogonal projection from the 
bundle TT> s \x>^ in TT> S to the tangent bundle TP* is given, at any 
rj G T>* and W G T V V S of the form W = w o rj for an H s vector field w 
on M, by the formula 

P(W) V = P(W o ?7~ 1 )id o 77, where P(u>)id = w - VA~ X divw. 

The fact that P is continuous on Ti&D s follows from the Hodge de- 
composition, but the fact that P is smooth as a function of the base 
point 77 is more subtle and requires the use of infinite-dimensional fiber 
bundles |EM] (we will use a similar technique to prove smoothness of 
the projection from TT> S E \jya to TT> s q in Theorem 13.41) . The resulting 
geodesic equation on "D s is the Euler equation 

D drj / drj dr]\ 
dt~dl ~ \dt'~dt)' 

written in terms of the covariant derivative of T> s and the second fun- 
damental form B; the more familiar form comes from writing ^ = wo 77 
and computing the orthogonal projection explicitly to obtain 

du 

— + V u u = -Vp, 

the Euler equation of ideal fluid mechanics. Smoothness of the con- 
nection implies that this is a smooth ODE on T"D* and leads to local 
well-posedness of the Euler equations from the Lagrangian point of 
view [EM] . 

We shall first study the induced connection on T> S E , the set of H s 
diffeomorphisms of M preserving both the volume form and a vector 
field E generating a circle action. By Corollary \2A\ this is a smooth 
sub manifold of T> s . If the field E is a Killing field for the given met- 
ric, the induced connection makes T> S E a totally geodesic Riemannian 
sub manifold of T> s . This implies that if a velocity field solves the Euler 
equation with initial condition commuting with E, then it will always 
commute with E. This is well-known in special cases such as axisym- 
metric Euler flow where E is an infinitesimal rotation about the axis. 
Of course, T> S E is not a totally geodesic submanifold of T> S E for the 
same reason D* is not totally geodesic in V s ; the only fluid flows that 
are also geodesies in V s are those for which the pressure function p is 
a constant function of x for each t, so Vp = 0. 
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Theorem 3.1. Suppose M is a compact Riemannian manifold with 
Killing field E with all orbits closed and of the same period. Then 
in the metric induced by ( ITT)) . the submanifold T> S E ^{M) defined in 
Theorem \2.3\ is a totally geodesic Riemannian submanifold ofD^(M). 

Proof. To prove this, we must show that whenever U and V are vector 
fields on T> E „, the covariant derivative VuV (computed by extending U 
and V to vector fields on D s and using the connection on TT> S ^) also lies 
in TT> S E . To make this easier, we convert from covariant derivatives 
of vector fields to covariant derivatives along curves: for any n G T> S E , 
let 7 be a curve such that 7(0) = 77 and '/(O) = U Tj . Define W to be 
the vector field along 7 given by Wit) = V^)- 

Now W(t) = w(t) o '-fit) for some time-dependent divergence-free 
vector field w on M, while U rj = u o rj for some divergence-free vector 
field u on M. By assumption [E, u] — and [E, w(t)] = for all t. The 
covariant derivative ^f- is given on by |Mi] 

Since [E 1 , = for all t, we can differentiate to get [E, ^|t=o] — 
as well, so we just need to show that [E, P(V u w)] = whenever 
divn = divw = and [E,u] = [E,w] = 0. 

Let £ be a diffeomorphism in the flow of E so ( : M — >• M is an 
isometry u] = u] = implies that C* M = u an d C*^ = ' u - 
Also since ( is an isometry, it preserves the covariant derivative. Thus 
C*V u i> = V u v. Furthermore, ( preserves the Hodge decomposition of 
vector fields, so £*PV u i> = P(*V u v = PV u v. Since this is true for 
any ( in the flow of E„ we find that [E, PV u v] = 0. Thus V s E fl (M) is 
totally geodesic in T> S (M). 

□ 

Let us now we use the connection induced on T> S E to define the 
connection on T> s q by orthogonal projection. As mentioned above, the 
hard part is showing that the orthogonal projection depends smoothly 
on the base point, but once we do this we get local existence of the 
geodesic equation for free. 

Lemma 3.2. Let M be a contact manifold of dimension In + 1 with 
contact form 9 and Reeb field E with closed orbits all of the same 
length, and suppose s > n + 3/2. There is a unique operator Sg from 
H S E + \M,R) to T ld V s E41 {M) such that for any f e H S E +1 (M,R), the 
vector field u = Sgf satisfies 9(u) = f and C u 6 = 0. 
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Furthermore, if Lg denotes the operator Lgu = C u 6, then Lg maps 
T id V E ^ into H S E ~ (A 1 ), the space of H s ~ l one-forms a such thata(E) = 
and LEda = 0. 

Finally, the sequence 

(19) H S E +1 (M,R) % T id V E ^M) H H E ~ 1 (A 1 ) 

is exact, i.e., the null space of each map is the image of the previous 
one. In particular the image of Sg is T i( iD s q {M) . 

Proof. We have already shown that So exists in Proposition II .61 To 
prove Sg has the desired smoothness, we use its expression in Darboux 
coordinates in (jjj). Since i?-invariance implies ^ = 0, we see Sgf is in 
H s if and only if / is in H s+1 . 

From above we have Lgu = d{9{uy\ + t u d9, so that if u is in H s then 
Lgu is only in fP -1 . The fact that Lg maps into ^^(A 1 ) follows as 
in the proof of Theorem 12.61 if [E, v] = and a = C v 6, then a(E) = 
and LEda = 0. □ 

Denote by Sg the formal adjoint of the operator Sg in the L 2 metric. 
To give an example, in Darboux coordinates with the Euclidean metric, 
it is easy to check that if w = (p k d x u + q k d y k) + r d z , then 

Formally, the orthogonal complement of the image of Sg is the null 
space of Sg, and the orthogonal decomposition of an arbitrary w G 
T\d^E^ is 

(21) W = PgW+(l-Pg)w, 

where Pg = SgAg Sq 
and Ag = SgSg. 

Again as an example we compute in Darboux coordinates with the 
Euclidean metric that the contact Laplacian Ag on a (not-necessarily 
^-invariant) function is 

(22) A 9/ . 2/ -£^(( 1 + ( ^)-£A|/. 

k=l k=l y y 

Observe that Ag differentiates twice in all directions except the z 
(Reeb) direction; this is a general feature. 

Note the similarity to the Hodge decomposition of a vector field into 
its gradient and divergence-free parts; this allows us to use the same 
methods to establish smoothness of the orthogonal projection as in 
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[EM] . There are now two issues: showing that the decomposition works 
for vector fields on M (which requires analyzing the smoothness of Sg 
and proving that Ag is invertible), then extending the decomposition 
from TidT>E,fi to TT>E tfl (which means by right-invariance that we have 
to prove that T R n PgT R~ l is smooth in rf). 

First we discuss the projection Pg at the identity. 

Theorem 3.3. Let M be a Riemannian manifold with contact form 
9 and a Reeb field E which is also a Killing field. The formal adjoint 
Sq of the operator Sg defined in Proposition \l.b\ and Lemma Iff.jjl maps 
T i( jP s E into H E ~ 1 (M,~R) , and the contact Laplacian Ag = SgSg is an 
operator mapping H E +1 (M,~El) bijectively onto H E ~ 1 (M,M). Hence the 
projection Pg defined by (I2TI) is a continuous map from TidD s E AM) to 
itself. 

Proof. We work one tangent space at a time: let x G M; we can then 
decompose T X M = RE © ker6» and T*M = W © ann£, where ann£ 
is the annihilator of E. Recall that the map 7 defined in Proposition 
11.11 is an isomorphism from ker# to axmE. Denote the inverse map 
from aimE to ker# by T, so that we have Sgf = f'E — Tdf for any 
/ G H E +1 (M,~R), since df annihilates E. Then the formula 

((S e f,w))= [ (fE-Tdf,w)dfi= [ f(E,w)-(df,rw)dfi 

JM JM 

= [ f(E,w) - fST*wdfx 

JM 

implies that SgW = (E, w)—5T*w, where 5 is the adjoint of d. Hence Sg 
is a first-order differential operator which maps T l( {D s into H S ~ 1 (M, M). 

To show that Sg maps T id V s E fl into H E ~ 1 (M, R), it is sufficient to 
show that if [E, w] = then E(Sgw) = 0. This will follow from the 
formula 

(23) SgCs = CeS$. 

Since £ is a Killing field, the operators Ce are skew-selfadjoint (on 
both functions and vector fields) which implies —CsSg = —S^Ce- The 
formula (|23|) is easily demonstrated (whether or not E is Killing) by 
writing Sgf in Darboux coordinates as in (JTj); clearly in coordinates Sg 
commutes with Ce = 

Finally to prove that* A e : #^ +1 (M,R) ->■ H S E ~ 1 (M 1 R) is invertible 
we recall that H E ±1 (M) is naturally identified with H S±1 (N) through 
f(n(x)) := f(x) where ir : M — > N is the Boothby-Wang fibration. 
With this identification we can consider A e : H S+1 (N) — > H^ 1 (N). 
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Furthermore we have 

(24) / fAgfdfi = [ \S e f\ 2 d^ = [ \fE-Tdf\ 2 d^ 
Jm Jm Jm 

for any ^-invariant function /. Also recall that for every x G M, we 
have T X M = RE © ker 9; using the fact that T is invertible and E is 
nowhere zero, we see that for each x there is a constant c x > such 
that 

(25) \f(x)E(x) - T x df(x)\ 2 > c x (f(x) 2 + \df(x)\ 2 ), 

and by compactness of M the numbers c x are uniformly bounded below. 

Since E is a Killing field on M, there is a unique Riemannian metric 
on iV such that the projection n is a Riemannian submersion. Hence 
the operator Ag is an elliptic operator on functions on N, and the 
formulas 021]) and fl2S) yield 

/ fA e fdv>c [ f+\df\ 2 dv, 

J N J N 

for any function / on N. (Here v is the induced volume form on N.) 
Thus identifying H^(M,M.) with H S (N,~R), we get a positive definite 
operator A e from H S+1 (N, M) to H S ~ 1 (N, R), which leads to an invert- 
ible operator from H S E +1 (M, R) to H S E \M, R) . □ 

Note that the contact Laplacian is never an elliptic operator on M, 
since it involves two derivatives in all directions except the Reeb di- 
rection; see (122]) . For example with the contact form on S 3 described 
in Example 11.3] we can easily compute that Sgf = fE\ — ^(E 3 f)E 2 + 
^(E2f)Es, and in the Berger metric we have 

(26) A e = a 2 - \E\ - \El 

On the space of functions with zero Reeb derivative, this causes no 
problem, but on the space of all functions on M, it requires dealing 
with a sub-Laplacian rather than a Laplacian. See [EbPrj for details 
on these issues. 

Having established that the projection makes sense at the iden- 
tity, we now want to extend it to a projection from the entire bundle 
TT> s E J\j)s (the set of all tangent spaces of T> S E ^ with base points in X>*) 
to TT> s q . Since the L 2 metric on T> E is right-invariant, the orthogo- 
nal projection is right-invariant as well: that is, for any rj G T> s the 

projection (Pg) v - T v T> s E fl -> T v V s q is given for any W G T V V S E/1 by 



(27) 



(PeUW) = (Pe(Wo V - 1 ))o V , 
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which is equivalent to (Pg) v = TR n o Pg o TR v -i . For this to be smooth 
on TT> S E M |x>| we need to establish that it is smooth in rj. 

Using the formula Pg = SgA^Sg, it is tempting to write 

(P e ) v = (TR v S e TR v -i) o (TR v Ag 1 TR r] -i) o (TR n S* e TR v -i) , 

expecting the three operators to be smooth functions of rj. In fact, 
for any first-order differential operator D the map TR V o D o TR v -i : 
H s — > H s ~ x is smooth in 77 by the Chain Rule; thus (TR v SgTR v -i) 
is smooth in 77. However the middle operator breaks down; for the 
whole projection to be continuous from H s to itself, we need the middle 
operator F 1— > (A^" 1 (F o r/ -1 )) o n to map if s_1 to H s+1 ; however the 
composition with 77 G T> s means that the result is only in H s , which 
breaks the entire result. Instead, we take an indirect approach following 
[EMj . Appendix A. We split the projection as 

(P e ) v = (TR n S e A^TR v -i) o (TR v S* e TR v -i) , 

and use the fact that the second map is smooth while the first is the 
inverse of the second on a restricted bundle (thus we obtain smoothness 
indirectly by the inverse function theorem). 

Theorem 3.4. The projection Pg given by ( J27l) is a C°° bundle map 
from the subbundle TT> S E to TT> s q . 

Proof. Since Sg is a first-order differential operator, the composition 
(JSf) v — TR V o SgoTR v -i is given on W G T v T> s Ell in a chart (using the 
chain rule) as a product of the inverse matrix of Drj and first derivatives 
of W . Since 77 is a diffeomorphism, the matrix Dr] is always invertible, 
and since rj is volume-preserving, Dr] has determinant one. Thus the 
map from Drj to its inverse is smooth, since by Cramer's rule it involves 
only determinants of minors, and products of H s ~ l functions are still 
in H s ~~ 1 because s > dimM/2 + 1 [EJ. We thus find Sg is smooth 
in any chart; hence it is smooth as a bundle map from TT> S E M (M) to 
V S E ^M) x H'-\M,R). 

Now consider TV s n as a subset of TV s Eii \ds. This is a subbundle 
since it is the null space of the bundle map Lg defined by (Lg) v = 
TR V o Lg o TR v -i as in Lemma T3.21 and Lg is smooth since as before it 
comes from a first-order differential operator. 

Since TV s q is a smooth subbundle, the restriction of Sg to TV S is 
still smooth. On this subbundle it is a bundle isomorphism onto its 
image V s q xH s E 1 {M, R), since TV s q is the image under the bundle map 
Sg of the bundle T> s q x H E +1 (M, M). Its inverse is given at each 77 by 
= TR n (S0Ag 1 )TR n -i. But the inverse of a smooth isomorphism 
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is also smooth by the inverse function theorem. Thus the map (Sg) 1 
is a smooth bundle map from V s q x H S E ~ 1 (M,R) to TV\. Since TV s q is 
a smooth subbundle, this bundle map is also smooth when considered 
as a map from V q x H s E l (M, R) into TT> S E fl \v a - Hence the composition 
Pg is also smooth as a bundle map. □ 

Now the orthogonal projection defines the Levi-Civita covariant de- 
rivative on the submanifold T> s q of T> E , and Theorem 13.41 implies that 
the induced Levi-Civita connection is C°°. This will allow us to discuss 
the geodesic equation in the next section. 

4. Local and global existence for the geodesic equation 

Since T> q is a Riemannian submanifold of T> S E , the covariant de- 
rivative on T> q comes from the orthogonal projection of the covariant 
derivative on T>% „, which we know is smooth. Thus the geodesic equa- 
tion on T> q is a C°° ordinary differential equation, which is guaranteed 
to have local-in-time solutions by the Picard method. 

Theorem 4.1. Suppose M is a Riemannian contact manifold, with 
contact form 9 and regular Reeb field E, such that E is Killing. Then 
for s > dim(M)/2 + 1 , the geodesic equation (P~e) (§^r) = on V q 
is equivalent to 

(28) ^L = S e fo V , ^Z + {/ ,A,/} = 0, 

where the contact bracket is defined for functions f and g in H E (M,M.) 
as in (JBJ) by {f,g} = S e f(g). 

The geodesic equation (|28|) is a smooth ODE on T> q , and hence there 
is a smooth exponential map exp id : Q C T id V q — > V s such that exp id (v) 
is the geodesic 7/(1) where Q is a neighborhood of zero in T id D q: 77(0) = 
id and r)'(0) = v. 

Proof. Let ^ denote the covariant derivative on T> E/X ; then a geodesic in 
T) s E tl satisfies = 0. Since V q is a smooth submanifold, its geodesic 
equation is (Pe) (£^) = 0- Having established smoothness of (Pe) 
in Theorem 13 A\ we know this geodesic equation is C°° on T> s . Since the 
geodesic equation is smooth, existence of an exponential map follows 
from the usual theory of ODEs in Banach manifolds. 

To compute the geodesic equation explicitly, we write ^ = u o r] 
for a time-dependent vector field u on M, which is in Ti&D q . We can 
write u = Sgf for some time- dependent function / in H E +1 (M). In 
the ambient space D s (M), we know the covariant derivative is given 
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^ dt ~dt = ( M * + ° an d the covariant derivative in V q is just the 
orthogonal projection onto TT> s q . To compute this, we let v = Sgg be 
an arbitrary vector in T- l( {D q , and define a right-invariant vector field 
by V n = v o 77. Then 

u) — (w, i>) divw — (w, [it, v}) + \\u\ 2 divi>) d/i. 

Now divw = divi; = since V q C I'b,^. Also [i/,i>] = [S0/, S^g] = 
Sg{f,g} where the contact bracket is given as in (jH]). Thus we find 

0= / (gA e f t - {/, g}A e f) d/i = [ g(A e f t + {/, A e f}) d/x. 

where the last equation comes from the standard property of Poisson 
bracket which was mentioned after the proof of Proposition 11.71 

Since this is true for any function g G i/^ +1 (M), we must have 

(HiE □ 

Example 4.2. Recall that for the contact structure on 5* 3 as in Exam- 
ple [L3l the contact Laplacian is given by (126]) : Ag = a 2 — \{E\ + E 2 ). 
As mentioned, this operator is more naturally understood in terms of 
the Boothby-Wang quotient, which in this case is iV = S 2 . We can 
then compute that the operator Ag reduces to a 2 — A on S 2 in terms 
of the usual Laplacian on S* 2 (the factor | cancels, since the metric on 
S 2 which makes the Hopf fibration a Riemannian submersion is \ of 
the usual one). Hence the geodesic equation (l2~8l) reduces to 

(29) <9,(A/-a 2 /) + {/,A/} = 

on S 2 in terms of the Poisson bracket (using the fact mentioned after 
Proposition 11.71 that the contact bracket on ^-invariant functions is 
equivalent to the symplectic bracket on the quotient). 

Equation (1291) is the quasigeostrophic equation on S 2 "in the /-plane 
approximation," where the parameter a 2 is the rotational Froude num- 
ber (or 1/a is the Rossby number) and / is the stream function; see 
Holm-Zeitlin |HZj and Zeitlin-Pasmanter [ZPj . (More precisely the f- 
plane approximation is valid only when the configuration manifold is 

^ote that the decomposition (|28p is valid is H s ~ 1 , but not in H s . The two 
equations are not actually ODEs on V s q , since left-translation is not differentiable 
in the H s spaces. The second formula in ([28]) is an expression of the general formula 
4| + ad* u = for the Euler- Arnold equation on a Lie group with right-invariant 
metric. This was originally derived by Arnold |Alj in the general case; see for 
example [A"K] or [MP] for more recent references. 
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IR 2 ; this is a simplification of the true quasigeostrophic equation.) In 
our context the Lagrangian motion is given in S 3 by ^ = Sgf o 77, so 
that 77 G V q (S 3 ). As we will see in a moment, the quantomorphism 
group is a circle bundle over the group of Hamiltonian symplectomor- 
phisms, and its Lagrangian flow projects to the usual Lagrangian flow 
on 5" 2 . We thus obtain a natural geodesic interpretation of the /-plane 
quasigeostrophic equation which is a bit different from the one given 
in jZP] . 

This works more generally; given any two-dimensional manifold iV 
with symplectic form u, we can find a three-dimensional contact man- 
ifold M such that iV is the Boothby-Wang quotient, and there is a 
unique choice of Riemannian metric on M such that the Reeb field E is 
orthogonal to ker 9 and |£^| 2 is a given constant Froude number; the ge- 
odesic equation on T> q (M) will then be the "/-plane" quasigeostrophic 
equation. Later (see Remark 14.81) we wills discuss the actual quasi- 
geostrophic equation on the 2-sphere. 

From the theory of ODEs on Banach spaces (see [L], Chapter 4, 
Proposition 2.5), we know that the only way a solution can blow up 
is if the norm blows up. Hence we can analyze global existence by 
establishing estimates for the Sobolev norms of u = Sgf in T ic {Dg(M). 
For this purpose it is convenient to view f[2"51) as an equation not on the 
contact manifold M but on its Boothby-Wang symplectic quotient N. 
Hence we consider Hamiltonian diffeomorphisms on N. 

Definition 4.3. We say that a diffeomorphism ( is in T>^ am (N) if 
there exists a curve ((t) from id to ( such that ((t) = u{t) o ((t) where 
u{t) = (*)) for some f(t) G H s+l (N) for every t. (Recall that u$ 
is the inverse of the operator : =04 i v w.) 

There is a natural fibration n from quantomorphisms T> S (M) on 
M to Hamiltonian diffeomorphisms 'Dl {aia (N) on N, first constructed 
by Ratiu and Schmid [RSj . The map is constructed as follows: if 
77 G T>g(M), then r\ preserves orbits of E by Proposition 11.21 Hence 
one gets ( G V S (N) by letting ((x) = tt(t](x)) where x = ir(x). Since 
r) G T> S JM), there exists a curve rj(t) from id to r\ in T> S {M) such that 
fj — vor) where v = Sgf for some curve / G H^~ 1 (M). Thus we can set 
u = n^Sgf = ir*(fE — Ydf) = — 7r*(rd/), and the flow £ of the time- 
dependent u will be a curve in T> S (N). But for any vector z G TM, 
we have u)(7r*rdf,7r*z) = d9(Tdf,z) = r yYdf{z) = df(z). Therefore 
u b {7iSdf) = df, where / G H S E +1 {M) is identified with / G H S+1 (N). 
Hence w b (vr,(?j)) = -df. It follows that ( G D^ am (iV). 
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Theorem 4.4. Let M be a Riemannian manifold of dimension 2n + 1 
with a regular contact form 9 for which the Reeb field E is a Killing field 
with orbits of constant length 1. Suppose that the Riemannian volume 
form fi on M is a constant multiple of A (d0) n . Define ip = A#(l) = 
Sq(E); then p is E -invariant and descends to a function (p: N — > IR 
such that <p o 7r = ip. Let N be the Boothby-Wang quotient with smooth 
projection tt: M — > N, and u the symplectic form on N satisfying 
n*u = d6. Then there is a unique Riemannian metric on N such that 
the projection it is a Riemannian submersion, and the induced volume 
form v on N satisfies n*v = and is a constant multiple of u n . 

Define a right-invariant Riemannian metric on T> s Ham (N) at the iden- 
tity by the formula 



(30) ((u,v)) N = / fAggdu p— , 

for vector fields u, v G T id V s Ham (N) , where u = u^df and v = u^dg, and 
Ag is identified with an operator on N as in the proof of Theorem \3.3[ 
Then the metric is well-defined on T> s Ham (N), and the map II is also a 
Riemannian submersion. 

Proof. Near any point of M we can choose a Darboux coordinate chart 
(x 1 , • ■ ■ , x n , y 1 , ■ ■ ■ , y n , z) such that 6 = dz + Y^k=i %k dy ] then the 
corresponding chart on N is (it 1 , ■ ■ • , u n , v 1 , ■ • • , v n ), and the projection 
in coordinates is 

vr(x\ • • • , x n , y 1 , ■ • • , y n , z) = (x 1 , • • • , x n , y 1 , ■ ■ • , y n ), 

i.e., just "forgetting about z." In such coordinates the Reeb field is 
E=-§- z , and since the Reeb field is Killing, the components of the Rie- 
mannian metric on M are independent of z. For each k the horizontal 
vectors are + p k and + q k ^ for some functions p k and q k 
(chosen to make the vectors orthogonal to E). Then if tt : M — > N is 
a Riemannian submersion, the metric on N satisfies 

(d uj , d u h) N = (d x3 + pP d z , d x k + p k d z ) M , 

(d uJ ,d v k) N = (d xJ + p 3 d z , d y k + q k d z ) M , 

(d vj ,d v k) N = (d y i + q 3 d z , d y k + q k d z ) M - 

The fact that the right side is independent of z implies that the left 
side is well-defined on N. 

Because both Riemannian metrics are right-invariant, it is sufficient 
to check the submersion condition at the identity. The construction of 
II shows that the null-space of II*: T^V^M) — > Ti^D^^N) consists 
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of vector fields Sgf such that df = 0. But for any such /, Sgf = fE 
with / constant, so the null-space is one-dimensional and spanned by 
E. Hence a vector Sgf G T id V q is horizontal if and only if ((Sgf, E)) = 
I'm / A e(!) dfi = J M fipdfi = 0. Now given any u G T id V^ am (N) there 
is a unique / : N — > E with lift / : M — > R satisfying f M fipdfi = 
such that u = oj*df; then U^(Sgf) = u, and so II is a Riemannian 
submersion if and only if ((u, u)) n = ((Sgf, Sgf))M for any such u. Since 
the orbits of E are assumed to have length 1, for any Reeb-invariant 
function if) on M we have J M ipdfi = f N ip dv for the induced function 
if) on N. 

To prove that ( l30l) is well-defined, we consider any vectors u and v 
in T id P^ am (A r ). Then u = u^df and v = u^dg for some functions / 
and g on N, which are defined only up to a constant. Now let / and g 
be the lifts of / and g to M, and write / = / + a and g = g + b, where 
we will demand that Sgf and Sg~g are horizontal. Then we must have 
f M fip dfi = f M g<p dfi = 0, which implies that 

a = -(Im f<P M/Um V and b = ~Um 9V dfi)/ (J M <p dfi). 

We compute that 

((Sgf, Sgg)) = J M fA g dfi + a J M gS e (l) dfi 

+ b f M fS e (l) dfi + ab f M Sgl dfi 
= Im f A e9 dfi - (J M fip dfi)(J M gip dfi)/(J M <p dfi) 
= In f A e9 dv - (J N fip du)(f N gip dv)/(f N ip dv). 

Hence although / and g are only defined up to constants, the metric 
fl30|) gives the same result for any choice of the constant and is thus 
well-defined as an inner product on the Hamiltonian vector fields u and 
vonV^jN). □ 

The metric (|30|) never reduces to the L 2 metric studied in [E2j , since 
that metric depends only on derivatives of the stream function. fl30|) is 
essentially a hybrid between the metric in |E2j and the L 2 Hofer metric 
on stream functions (which only generates a pseudometric on the space 
of Hamiltonian diffeomorphisms; see |ElPoj .) In the situation described 
in Example 14. 2| the induced metric (1301) on symplectic diffeomorphisms 
of S 2 is exactly the H 1 metric on stream functions. We must be careful, 
however: since the tangent space of V^^N) consists of skew gradients 
of functions, it can only be identified with functions defined up to 
some constant. Hence the metric (!30l) is necessarily degenerate in one 
direction. Nonetheless we can still write a geodesic equation as an 
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Euler-Arnold equation on the corresponding homogeneous space (see 
Khesin-Misiolek [KM] ). 

Proposition 4.5. On the space T> s Ham (N) with right-invariant degen- 
erate Riemannian metric given by fl30|) . the Euler-Arnold equation on 
the homogeneous quotient space is given by 

(31) d t A e f + {f,A e f} = 0, 

where {■,■} is the symplectic Poisson bracket and f is determined only 
up to addition by a constant. 

Proof. We just need to compute the Euler-Arnold equation u t +ad* u = 
for u = urdf in the metric ( J30l) . Using a.d u v = —u^{f,g} for u = 
ufidf and v = oo^dg, we can easily compute that the Euler-Arnold 
equation is 

9Aef ~ S ~i^ + {f, A e f} - ^{f, ^ = 0. 

Choosing a representative / such that j N f(pdv = 0, we obtain ( 13T|) . 

□ 

The geodesic equation fTSSl) is equivalent to the equation (jHTj) for a 
time-dependent function / : N — > R, where is the inherited operator 
defined as in Theorem 14.41 Global existence for the quantomorphism 
geodesic equation then reduces to showing that this equation has global 
solutions for sufficiently smooth initial data: since we want Sgf to ini- 
tially be in H s , we assume that / G H s+l (M, R). Global existence 
follows in much the same way as in [E2j : we first bound the Sobolev 
norms of / in terms of the C l norm, then establish an absolute maxi- 
mum for the C 1 norm. In that paper the proof was very brief; here we 
give more detail. 

Theorem 4.6. If the initial condition f is in H^ 1 (M, R) for s > 
dimiV/2 + 1 (i.e., if f is in H S+1 (N,'R)), then solutions of equation 
(|3T|) exist for all time. Hence the exponential map given by Theorem 



4-l\ is defined on the entire tangent space TidD^(M). Thus by right 
invariance it is defined on all ofTT> s q (M). 

Proof. In order to prove this theorem we recall that the space if^, ±1 (M, R) 
is naturally identified with H S±1 (N, R). With this identification we get 

A e : H S+1 (N,R) -> iT-^R) 

We shall deal with three time dependent quantities whose estimates 
will be interdependent. We have the stream function f(t) G ifg +1 (M, R) = 



28 



DAVID G. EBIN AND STEPHEN C. PRESTON 



# s+1 (iV,M), the velocity field u(t) = S f(t) e T id V s q (M) and the La- 
grangian flow rj(t) in T> S JM). Let (— T e ) be the maximum time 
interval on which / (and therefore u and 77) exists. If T e were finite, 
then we would have \\m t/ ^T e ||/(£)|| s +i = 00, because if ||/(£)|| s+ i re- 
mained bounded then ||w(t)|| s would be bounded and then T e would 
not be maximal (see [H]). Hence we need only show that is 
bounded on any finite time interval. 
To do this we first note that 

^-(A e f(t) o v (t)) = (d t A e f + S e f(A e f)) o V = 
at 

by d2HD. Therefore A e f{t) o rj(t) = A e f(0), so in particular the C°- 
norm of Agf(t) is constant in time. Using this we will find a bound for 
df, or equivalently for Sgf. 

Recall from the proof of Theorem 13.31 that the contact Laplacian 
Ag: H S+1 (N,M.) — > H S ~ 1 (N,M.) is an isomorphism. From the theory 
of elliptic operators (PQ, Chapter 7, Proposition 2.2 for N = MJ 1 and 
Chapter 7, Section 10 for N any compact manifold) we find that its 
inverse is a pseudodifferential operator whose Schwartz kernel k : N x 
N — > K is smooth off of the diagonal and obeys the estimate 

(32) \d£k(x,y)\ < Kp{x,y)- 2n+2 - M for - 2n + 2 - \(5\ < 

where d x means a partial derivative with respect to the first variables 
of k and p(x,y) is the distance from x to y defined by the Riemannian 
metric on Nu From this if n > 1, we can estimate 

f(x)= [ k(x,y)A e f(y)dy 

J N 

and 

(33) df(x)= / d x k(x,y)A e f(y)dy 

J N 

where d x is the differential with respect to the ^-variables and dy in- 
dicates integration with respect to y using the Riemannian volume 
element of iV"|l We have \k(x,y)\ < K p(x,y)~ 2n+2 and \d x k(x,y)\ < 
K p(x , y)~ 2n+1 , so these integrals are bounded by a constant times 
|| Ag/^t) llc-o which, as we have seen, is independent of t. With this 
and formula ([7]) we see that u = Sgf is bounded uniformly in time 



In the inequalities that follow K will always denote some positive constant, but 
it may different in different inequalities. 

3 In the sequel we will sometimes write dk for d x k. 
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as well. If n = 1 the estimate for df is still valid, but we must use a 
different method to estimate /. First we note that if \df | < K then 

(34) max / — min / < KR 

where R is the diameter of N. Also we see from (I3TI) that J N Agf is 
constant in time. Furthermore 

/ Agf = [ fA e l = [ ftp 

JN JN JN 

and I N f = Iai\ E \ 2 > °- Also Inff ^ fA min f)¥ = min / Jn<P- 
Thus min / is bounded so by fl34|) / is also bounded, and of course all 

the bounds are uniform in time. 

We proceed to seek a time-uniform Lipschitz bound for df, but in 
fact we will be able to find only a quasi-Lipchitz bound, as we shall 
now explain. Since A" is compact, we can find a positive 5 such that 
each point of A" has a normal coordinate neighbourhood ball of radius 
at least S. Then for y G A" with p{x, y) < 5, we have a unique minimal 
geodesic x parametrized so that x(0) = x and xiX) = V- We let 
b = |x'(r)|, so that b = p(x,y). Also we parallel translate df(x) along 
X to get some df'(y) G T*N. Then we shall estimate \df(y) — df'(y)\ 
where | | is the norm on T*N. 

To do this we use the formula (1331) for df, and we parallel translate 
each d x k(x, z) along x to get dk'(y, z) G T*N. In this way we get 

df'(y)= / dk'(y,z)A e f(z)dz 

JN 

and we find 

df'(y) - df(y) = / (dk'(y,z) - dk(y,z))A e f(z)dz 

JN 

Following [Kaj . Lemma 1.4, we divide this integral as follows: Let 
£ = B 2 b(x), the ball of radius 2b about x. Then f N = J E + For 
the integral over S, we have the estimate 

/ dk'(y,z)A e f(z)dz<K\\A e f\\ c o [ p(x, zf- 2n dz < 2bK\\ A e f\\ c o. 

JT, JT 

Also E C B 3 b(y) and 
/ dk(y,z)Agf(z)dz<K\\A e f\\co [ p(y, z) l ~ 2n dz < 3bK\\ A 6 f\\ c o 

JT JT 

Thus 

(35) / (dk'(y, z) - dk(y, z))Agf{z) dz < 5bK\\A e f\\ c o 

JT 



30 DAVID G. EBIN AND STEPHEN C. PRESTON 

The estimate of the integral over N — £ is more subtle. If we let 
P(t) denote parallel translation along x from xi T ) to x(fy, so that 
P(r) : T* (t) N -> T y *A^, we find that 

(36) dk'(y,z) -dk(y,z) = P(b)(dk(x, z)) - dk(y, z) 

= [ P(T)V x/{T) dk( X (r),z)dz 
Jo 

Therefore 

(37) \dk'(y,z)-dk(y,z)\ < max{| V x , {T) dk{x{r), z)\} 

0<t<o 

< Kb{p{x,z)-b)- 2n 
We note that R = diam(iV). Hence 

(38) f \dk'{y,z)-dk{y,z)\ < Kb! (p(x, z) - b)~ 2n dz 
Jn-t, Jn-y, 

f R 

< Kb (r-by 1 dr 

J2b 

= Kb(log(R-b)-\ogb) 

< Kblog(Rfb) 

Combining this with (1351) we find 



A' 



\dk'(y, z) - dk(y, z)\dz< K(5p{x, y) + p(x, y) log(R/p(x, y)) 
for p(x, y) < 5. By increasing K we simplify this inequality to 
\dk'(y,z) - dk(y,z)\dz < Kp(x,y)(l + \og(R/p(x,y)). 



N 



Thus we find that if p(x, y) < 5, 

(39) \df({y)-df(y)\<Kp{l + log(R/p(x,y)) 

where K is independent of t as before. With this inequality we say that 
df is quasi-Lipschitz. Also since df(x) is bounded independently of x 
and t we find that by further increasing K we get ( 139]) for all x,y G N; 
that is, we can drop the restriction p(x, y) < 5. 

Since u = Sgf, our bound for df gives the same quasi-Lipschitz 
bound for u, again uniformly in t. Using this bound we can find a 
positive a for which the flow rj(t) of u is C a . However we will need 
to show that rjit)^ 1 is C a . Fortunately we can do this by the same 
method. 

Given t £ [0, T e ) we define a time dependent vector field v on M 
by v(t) = —u(to — t). Then if a is the flow of v, it is easy to see 
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that the maps a(t, rj(to, x)) and 77(^0 — t, x) satisfy the same differential 
equation, and since a(0,r](t ,x)) = r](to,x), they must be equal for all 
times t G [0, to]- Hence in particular a (to) = 77(to) _1 - We proceed to 
show that a (to) is C a . 

Fix x and y in M. Let x : [0,to] x [0, 1] be the map such that r 1— >■ 
x(t, t) is the minimal geodesic between a(t, x) and a(t, y) with x(t, 0) = 
a(t,x) and x(t, 1) = cr(t,y). Define 



(f)(t) = p(a(t,x),a(t,y)) 

Then 



./o 



*m = £(/W>i*- 



1 



'MM,^(t,r))dr 



/o |gf(t,r)| 

But J^l^ = by general properties of surface maps, and = 
since each r !->■ x(t, t) is a geodesic, and thus an integration by parts 
yields 



1 



|T=1 



using the fact that 1 1^ (t, r) | is constant in r since x is a geodesic in r. 

Now |^ is parallel along x, and since the parallel transport P from 
x(0) to preserves inner products, we have 

<§* (t, 0), a(t, x))) = (§* (t, 1), a(t, a:))). 

Thus 

< -Pv(t,a(t,x))\. 

Now since w like u is quasi-Lipschitz on [0, to], we have a constant .ff 
such that 

(40) 0'(t)<^(t)(l + log^^, 

where the constants R and K do not depend on t or t . 

We proceed to estimate <f)(t). Let ^(t) = log(<f)(t) / R) so = 
Then from (HO)) we get ^' < — ?/>). Integrating this we find 

ip(t) < ij(0)e~ Kt + 1 - e~ Kt . 
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Exponentiating this inequality and noting that (f)(0) < R we find that 

m < rm \ Kt e i-.-« < (m 

R ~ V R ) ~\R 

Thus <j)(t) < R 1 - e ~ KTe e(j)(0) e ~ KTe , so letting a = e~ KT - and L = 
eR}~ e e , we get <f>(t) < L<f)(0) a and hence 

(41) p(a(t , x), a(t , y)) < Lp(x, y) a . 

The constants L and a do not depend on the choice of to, so the estimate 
fj4~lj) holds for all t <G [0,T e ). From (JH]) and the equation 

A fl /(t) = A fl /(0)o<r(t) 

we find that Agf(t) is bounded uniformly in t in C a (M, R), and thus 
by standard elliptic theory we get f(t) bounded in C 2+a (M,~R), from 
which it follows that Sgf is bounded in C 1+a (TM). 

Now we show f{t) to be bounded in the H s+1 topology, or equiva- 
lently that Agf(t) is bounded in H s ~ l . We note that 

d t A e f = -S e f(A e f) 

so 

i / {A e ffp = - ! Sof(A e f) 2 p= [ {A e ff div(S e f)p = 
Jm Jm Jm 

Taking s — 1 spatial derivative^ we get 

/ \V s - 1 A e f\ 2 p = -2 [ (V 5 ./V- 1 A«/,V- 1 A tf /)A* 

-2 / ([V5,/,V- 1 ]A fl / ) V- 1 A tf /> A i J 

where [ , ] denotes the commutator. The first term on the right van- 
ishes, and for the second we use the standard estimate 

\\V l (hg) - hV l g\\ H o < K (\\h\\ H i\\g\\ c o + ||V/i|| c o||#||^-i), 

which can be found in [T] Chapter 13, Proposition 3.7. We let I — s — 1, 
h = Sgf and g = Agf and obtain 

/ (\ysof, V s ~ l ]A e f, V s ~ l A e f)p < K\\Sof\\c4*of\\%.-i. 
Jm 

4 Powers of V are defined in a standard way: for any function /, V/ is a section 
of TM, V 2 / is a section of TM ® T*M and for any k, V fe J is a section of TM ® 
(T*M) fe_1 . Inner products are defined by the induced Riemannian metric. 




33 



Thus we get 



j t \\^ef\\ 2 H s-i < K\\Sef\\ c 4A e f\\ 2 Ha 



from which it follows that ||Ae/||^-i remains bounded on [0,T e ), since 
ll^/Hc is bounded uniformly in time. The same reasoning combined 
with a time reversal gives boundedness on (— T&, 0]. This completes the 
proof. □ 

Remark 4.7. We also get C°° geodesies for all time. If the initial datum 
is C°°, the geodesic exists in V s q . But for any s' > s, there is also a geo- 
desic in D s ' C T> s . By uniqueness of solutions for ordinary differential 
equations, the two geodesies must coincide. Since s' can be arbitrarily 
large the geodesic must be C°°. 

We have shown in Example 14 . 2 1 1 hat the geodesic equation on T> q (M) 
reduces, with a certain choice of metric on M, to the equation 

(Af-a 2 f) t + {f,Af} = 1 

which is a simplified version of the quasigeostrophic equation; when 
M = R 2 , this is the /-plane approximation. Now we finally want to 
discuss the real quasigeostrophic equation on S 2 , which is 

(42) (A/ - a 2 f) t + {/, Af-fi cos 9} = 

for some parameters a and /3. Here we use spherical coordinates in 
which z = cos 9 is the height function on the 2-sphere. See |Kul ICS} 
ISTSj for derivations of this equation. The equation (1421) can be written 
in conservation form in terms of the "potential vorticity" u = Af — 

(43) u t + u{u) = 

where u = Sgf. Now let us show how to derive this equation in our 
geometric context as an Euler-Arnold equation. 

Remark 4.8. Consider the Lie algebra Ti&D q (M) = J 7 e(M), with bracket 
given by [Sgf,Sgg] = Sg{f,g}, where {■,■} is the contact bracket de- 
fined by ([Bj). Let tp: M — > R be an arbitrary smooth function, and 
construct a central extension structure on the product Ti&D q (M) x R 
by the rule 

(44) [(S e f,a),(S e g,b)} = (s 9 {f ig },j^{f,g}dy\ . 

It is easy to check that this bracket satisfies antisymmetry and the 
Jacobi identity, so that it is a genuine Lie bracket on a Lie algebra. 
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Now define an inner product on T\ ( {D q {M) x R by 
(45) (((S 6 f, a), (S e g, b))) = [ (S e f, S e g) d» + ab. 

J M 

Then it is straightforward to compute that the Euler- Arnold equation 
for this Lie algebra and inner product is the system 

A e f t + {f,A e f} + a{f,<p} = 0, 

a t = 0. 

In particular if M = S 3 with Berger metric as in Example (14. 2p . and 
if tp is the height function on the sphere S 2 , then we can obtain the 
genuine quasigeostrophic equation on S 2 with an appropriate choice of 
parameter a. 

The proof of global existence in Theorem 14. 61 can be extended without 
much difficulty, since everything there depends on having a uniform 
bound for the potential vorticity Agf — a 2 f, which now gets replaced 
by the real potential vorticity Agf — a 2 f — (3cos0. Obviously a uniform 
bound on one of these implies a uniform bound on the other, so that all 
the techniques immediately generalize, and we obtain global existence 
for solutions of the real quasigeostrophic equation on S 2 . 
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